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ABSTRACT 



Context. Observations with the Hinode spectro-polarimeter have revealed strong horizontal internetwork magnetic fields 
in the quiet solar photosphere. 

Aims. We aim at interpreting the observations by means of results from numerical simulations. 

Methods. Radiative MHD simulations of dynamo action by near-surface convection are analyzed with respect to the 
relation between vertical and horizontal magnetic field components. 

Results. The dynamo-generated fields show a clear dominance of the horizontal field in the height range where the spec- 
tral lines used for the observations are formed. The ratio between the averaged horizontal and vertical field components 
is consistent with the values derived from the observations. This behavior results from the intermittent nature of the 
dynamo field with polarity mixing on small scales in the surface layers. 

Conclusions. Our results provide further evidence that local near-surface dynamo action contributes significantly to the 
solar internetwork fields. 
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1. Introduction 



The ubiquitous existence of small-scale 'internetwork' 
magnetic fields of mixed polarity in the so-called 
quiet solar photosphere is strongly indicated by vari- 
ous observational diagnostics ( e.g., Khomcnko ct al.,, 
iLites fc Socas-Navarrol 12004 iTruiillo Bueno et all , 
and further references therein). Recent high- resolution 



2. Numerical model 



2003; 



2004L 



space-borne observations with the spectropolarimeter of 
the Solar Optical Telescope aboard the Hinode satel- 
lite have considerably strengthened the case for inter- 
network fields and, furthermore, have revealed that the 
measured internetwork flux is dominated by strongly 
inclined, almost horizontal mag netic field (|Lites et al.l . 
120071 10 rozco Suarez et all . l2007bh . Considerable amounts 
of highly time-dependent horizontal magnetic flux have also 
been found in g round-based observ ations with lower spa- 
tial resolution ((Harvey et al.l . l2007t ). The ubiquity of the 
small-scale, mixed-polarity internetwork field suggests a lo- 
cal origin of at least a significant part of the measured 
flux. Recently, we have demonstrated by means of radia- 
tive magneto-convection simulations that local dynamo ac- 
tion by near-surface co nvective flows is a possible source for 
the internetwork flux (|Vogler fc Schiisslerl |2007[ ) . Here we 
show that the spatial structure of the dynamo-generated 
field provides a natural explanation for the observed domi- 
nance of the horizontal field component in the middle pho- 
tosphere. 



We us e the results of dynamo run C of IVogler fc Schiisslerl 
()2007D with 648 x 648 x 140 grid cells in a computational 
box with a physical size of 4.86 x 4.86 Mm^ in the horizon- 
tal and 1.4 Mm in the vertical direction, the latter ranging 
from about 900 km below to 500 km above the average 
level of continuum optical depth unity at 630 nm wave- 
length (rgao = 1). The simu l ation has been run with the 
MURaM code (|Vogled . [20031: IVogler et al.l . [200l . With a 
magnetic Reynolds number of about 2600, the simulation 
shows an exponential growth of a weak seed field with an 
e-folding time of about 10 minutes. The magnetic energy 
saturates at about 2.5% of the kinetic energy of the convec- 
tive flows, the maximum of the spectral energy distribution 
lying at horizontal spatial scales of a few hundred km, at 
which scales the field displays a distinctly mixed-polarity 
character. 

3. Relation of horizontal and vertical field 
components 

Figure [T] is based upon a snapshot from the saturated phase 
of the dynamo run. The intensity image in the upper left 
panel shows the granulation pattern, which is almost undis- 
turbed by the presence of the magnetic field. The distri- 
bution of the vertical field component on the level sur- 
face Tgao = 10~^ (upper middle panel) reveals the mixed- 
polarity nature of the dynamo-generated magnetic field, 
which preferentially resides in the intergranular downflow 
lanes. Note that the distribution of vertical field at this level 
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Fig. 1. Snapshot from the saturated phase of the dynamo run C (j Vogler fc Schiissleij . l2007[) . The panels in the upper row 
show the continuum intensity at 630 nm wavelength (upper left), a gray-scale image ('magnetogram', saturated at ±10 G, 
black and white colour indicating the two polarities) of the (signed) vertical field component on the surface Tgso — 10^^, 
roughly corresponding to the formation height of the spectral line used for the observations (upper middle), and the 
unsigned horizontal field strength (upper right), black colour in this panel representing very weak field and white colour 
indicating the saturation level of 10 G. The two panels in the lower row show gray-scale images of the logarithm of the 
horizontal field strength (black colour indicating fields below 1 G) for vertical cuts through the simulation box along a 
horizontal line at y = 1.5 Mm (left) and y — 3 Mm (right), respectively, in the upper panels. The level of Tg3o — 1 is 
roughly at a height of 0.9 Mm, the level of Tgso = 10^^ at 1.2 Mm. 



is significantly smoother than at the height of optical depth 
unity (cf. Fig. 2 of .Vogler fc Schii ssler. 2007). This reflects 
the fact that much of the small-scale flux at the lower level 
has already been connected back by shallow loops. As a 
consequence, the unsigned vertical field drops much more 
rapidly with height than the (unsigned) horizontal field, 
so that at the level rgso = 10~^ the latter (shown in the 
upper right panel) dominates in most places. The represen- 
tation on the two vertical cuts shown in the lower row of 
Fig.[T]illustrates that the horizontal field in the photosphere 
(i.e., above a height of about 0.9 Mm) has two components, 
namely, narrow loops near the intergranular lanes and ex- 
tended loops above gr anules. The latter have also been seen 
in the s imulations of iGrossmann-Doerth et all (|1998 ') and 
[St eined (2007); they are presumably formed by reconncction 
events between loop 'legs' with opposite polarities and flux 
expulsion by the granular flows, larger (stronger) granules 
pushing the horizontal field to higher levels in the photo- 
sphere. As a consequence, a cut at a given level surface of 
constant height (or optical depth) as shown in the upper 
right panel of Fig.[l]misses part of the horizontal field above 
granules and will be m ore dominated by the field around 
the intergranular lanes. iLites et al] (|2007f) report that the 



observed horizontal field is spatially separated from the ver- 
tical field and favors the edges of bright granules. A com- 
parison with this finding requires an analysis based upon 
synthetic Stokes profiles from the simulation data, which 
will be presented in a later paper. 

A quantitative account of the relation between verti- 
cal and horizontal field components is given in Fig. [2l It 
shows field strengths averaged over surfaces of constant op- 
tical depth as functions of rg3o . The dashed curve gives the 
average unsigned vertical field while the dotted and dash- 
dotted curves represent the averages of the two horizontal 
field components (unsigned). 

The averages of the three components indicate that the 
field is not far from being statistically isotropic in the deep 
layers below Tgso = 1. In contrast, the horizontal compo- 
nents of the magnetic field become increasingly dominant 
in the photosphere above. The driving of the dynamo by 
small-scale turbulent shear flows in and adjacent to the in- 
tergranular downflows is mainly restricted to the regions 
below T630 = 1., whereas, in the convectively stable photo- 
sphere above, these flows are much weaker and the rate of 
work against the Lorentz force drops steeply with height. 
Since inductive effects have smaller influence on the field 
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Fig. 2. Profiles of magnetic field strength averaged over 
surfaces of constant rgso: unsigned vertical field (dS^I), 
dashed), unsigned horizontal field components {{\Bx,y\): 
dotted and dash-dotted, respectively), and root-mean- 
square of the horizontal field strength ((B^ + 5^)^/^, solid). 



Fig. 3. Ratio of the root-mean-square of the horizontal field 
component (solid line in Fig. [2]) to the averaged unsigned 
vertical field (dashed line in Fig. In the optical depth 
interval —2 < log(T63o) < — 1 (relevant for the formation of 
the Fel lines used in the Hinode SP), the ratio is roughly 
in the range 4-6. 



in the layers above Teao = 1, the decay of the field with 
height is mainly determined by its spatial structure at the 
surface (particularly by the energy spectrum as a function 
of horizontal wavenumber). This results in a steep decline 
of the unsigned vertical field with height as opposite polar- 
ities on small scales are connected by shallow loops with 
typical length scales of a few hundred km, corresponding 
to the horizontal scale for which the magnetic energy spec- 
trum at Te3o = 1 reaches its maximum. It becomes plausible 
that this configuration leads at the same time to a less steep 
decline of the horizontal field, if one considers the simple ex- 
ample of an arcade-like magnetic field with concentric semi- 
circular field lines: for increasing height, more and more 
field lines turn over horizontally, so that the horizontally 
averaged unsigned vertical field strength decreases faster 
than the averaged horizontal field; a simple calculation for 
this case shows that (|-Bhoi |) strongly exceeds (|-Bvcrt|) at 
heights of the order of the horizontal scale (footpoint sep- 
aration at the surface) of the arcade. Thus, the dominance 
of horizontal fields in the photosphere is consistent with 
the assumption of a simple loop topology with a preferred 
length scale. 

The quantity that is actually relevant for a qualitative 
compari son with the 'apparent' horizontal field strength de- 
rived bv lLites et al.l ([^007) from measurements of the linear 
polarization (transversal Zeeman effect. Stokes Q and U) 
is the root-mean-square of the horizontal magnetic field, 
i.e., -Brms = {Bl + By)^/"^, since for not too strong fields 
Stokes Q and U are proportional to the square of the hor- 
izontal field strength. Srms as a function of optical depth 
is shown as the solid line in Fig. [21 Owing to the inhomo- 
geneity of the horizontal field, this quantity is significantly 
larger than (I-BkI) and (|i?y|), even if we multiply any of 
them by a factor to take into account both horizontal 
field components. From the ratio of the average horizon- 
tal field and the rms field, we can estimate an average 'fill 
fraction' as a measure of the inhomogeneity of the hor- 



izontal field. We obtain a number between 0.25 and 0.3 
in the range —2 < logrgso < —1, which is the relevant 
range of line formation of the Fel lines at 630.15 nm and 
630.25 nm used by t he Hin ode spectro-polarimeter (e.g., 
lOrozco Suarez etall . l2007af) . These values ar e consistent 
with t he es timate of the fill fract i on obta ined bv lLites et al.l 
()2007D and lOrozco Suarez et~all (|2007bt) by means of an in- 
version method. 

Let us now consider the ratio between the average ver- 
tical field and -Brms as shown in Fig. [S] The ratio increases 
strongly with height, so that in the optical depth interval 
—2 < logTgao < — 1 values between 4 and 6 are reached. 
This is consistent with the ratio of the horizontally av- 
eraged vertical (longitud inal) and horizont al (transversal) 
apparent fields found bv iLites et al.l (|2007D : Sfpp/S^p = 
55 G/11 G = 5. However, as pointed out by these authors, 
the relation between i3?l„ and the actual horizontal field 
strength is far from trivial. Furthermore, effects of line-of- 
sight integration and spatial smearing by the instrument 
complicate the relationship between the average fields in 
the simulation and the field strengths derived from the 
observed Stokes profiles. A direct quantitative comparison 
with the observations would have to proceed by means of 
calculating synthetic Stokes profiles, taking into account 
the point-spread function of the instruments. This is be- 
yo nd the scope of thi s Letter. 

iLites et all ()2007D have suggested that one possibility 
contributing to the imbalance of the average vertical and 
horizontal fields could be a significantly larger horizon- 
tal scale of the horizontal field as compared to the verti- 
cal field. In fact, this is what we clearly find in our dy- 
namo simulation. Fig. [4] shows spectral magnetic energy as 
a function of horizontal wave number. The dashed curve 
gives the energy distribution for the vertical field, while 
the solid curve represents the spectral energy in the hor- 
izontal field (mean of the spectra for the two horizontal 
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Fig. 4. Magnetic energy spectra as a function of horizontal 
wave number, fc, for the height range roughly corresponding 
to — 2 < logrgso < — 1. The dashed curve shows the energy 
spectrum based on the vertical field component, while the 
solid curve gives the arithmetic mean of the spectra for the 
two horizontal field components. For low wavenumbers, the 
energy in the horizontal field clearly dominates. 

field components). For this plot we have considered fields 
in the height range roughly corresponding to the optical 
depth interval — 2 < logreao < — f, which is relevant for 
the formation of the iron lines used for the observations. 
The curves show that the field components are in equipar- 
tition at small scales (large wave numbers), but that the 
horizontal field clearly dominates at wave numbers below 
roughly fO Mm~^, corresponding to horizontal scales larger 
than about 600 km. 

4. Discussion 

A direct comparison of the simulation results with the ob- 
servations is not possible because the individual values for 
the averaged fields in the simulation are both (by about a 
factor three) smaller in the relevant height range than the 
values for the 'apparent' derived from the observation, the 
discrepancy probably becoming even more severe when the 
actual spatial resolution of the observations is taken into ac- 
count. This is not particularly surprising since the magnetic 
Reynolds number of the simulation is still orders of mag- 
nitude smaller than the actual value in the corresponding 
solar layers, so that the saturation level of our simulation is 
probably considerably below the level to be expected for the 
real Sun. In fact, a preliminary simulation with a roughly 
doubled Reynolds number of about 5000 shows an increase 
of the magnetic energy by a factor of about 1.7 with respect 
to the case shown here. Interestingly, it turns out that the 
optical depth profiles of the averaged field strengths in this 
case can be very well approximated by just multiplying the 
curves shown in Fig. [2] by Vl.7, meaning that the main dif- 
ference between the simulations reduces to a simple scale 
factor for the field strength. Together with the fact that 
the observed internetwork fields are predominantly weak 



and their energy is significantly smaller than the kinetic 
energy of the convective motions, this suggests that the 
general nature of the dynamo-generated field may not be 
significantly different from the case shown here, apart from 
a higher overall amplitude. In particular, we expect that 
the ratio of the average horizontal and vertical field is not 
strongly affected by the amplitude of the dynamo-generated 
field. Of course, these assertions need to be demonstrated 
by further simulations with higher Reynolds numbers. 

The clear dominance of the horizontal field in the mid 
photosphere seems to be a rather specific property of the 
strongly intermittent field generated by near-surface tur- 
bulent dynamo action . Accordingly, the dynamo simula- 
tion of lAbbettl (|2007l with a closed bottom boundary 
and a local treatment of radiative transfer) also exhibits 
strong horizontal field in the photospheric layers. On the 
other hand, models with an imposed net vertical flux (e.g., 
iVogler et all . [2005h or our recent simulations of the de- 
cay of a granulation-scale mixed-polarity field (at magnetic 
Reynolds numbers below the threshold of dynamo action) 
do not show this behavior; in these cases, the intricate 
small-scale mixing of polarities that is characteristic for the 
dynamo does not dominate the field structure. The rapid 
decay with height of such a dynamo-generated field is also 
consistent with the a pparent lack o f stron g horizontal field 
in the chromosphere (|Harvev et al.l . [20071 ). 

What are the alternatives to near-surface dynamo ac- 
tion? 'Shredding' of pre-existing magnetic flux (remnants of 
bipolar magnetic regions) cannot explain the large amount 
of observed horizontal flux since the turbulent cascade does 
not lead to an accumulation of energy (and generation of 
a spectral maximum) at small scales. On the other hand, 
such a behavior is typical for turbulent dynamo action. Flux 
emergence from the deeper convection zone in the form of 
granule-sized small bipoles would have to proceed such a 
high rate in order to maintain the ubiquitous strong hor- 
izontal fields that it probably would not have gone un- 
detected in the past (see, however, iCenteno et al.l . 120071 ). 
The sporadic appea rance of hor i zonta l internetwork fields 
(HIFs) described by iLites et all (|1996f ) and interpreted as 
small-scale flux emergence events seems to be unsufhcient 
to explain the ubiquitous horizontal field now found with 
Hinode. On the other hand, flux recycling of an overall back- 
ground flux by granulation probably represents a signifi- 
cant source of horizontal field in network and plage regions. 
Emergence of e xtended horizon t al fiel d strands in granules 
as observed by llshikawa et al.l (|2007f ) is not seen in local 
dynamo simulations. 

In the real Sun, probably all three sources, i.e., dy- 
namo, shredded fields, and small-scale flux emergence from 
deeper layers, contribute to the internetwork flux in un- 
known amounts. In any case, the strong horizontal fields 
in the quiet photosphere inferred by the observations in- 
dicate that the source of these fields at the solar surface 
is a mixed-polarity field whose energy is mostly contained 
in those spatial scales where the dynamo-generated flux 
resides. Therefore, the observational results obtained with 
the Hinode SP together with the analysis presented here 
provides strong indication that surface dynamo action rep- 
resents a significant source for the internetwork field in the 
solar photosphere. 
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